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Abstract. We employ an adaptation of a strong-disorder renormalization-group 
technique in order to analyze the ferro-paramagnetic quantum phase transition of Ising 
chains with aperiodic but deterministic couplings under the action of a transverse field. 
In the presence of marginal or relevant geometric fluctuations induced by aperiodicity, 
for which the critical behavior is expected to depart from the Onsager universality 
class, we derive analytical and asymptotically exact expressions for various critical 
exponents (including the correlation-length and the magnetization exponents, which 
are not easily obtainable by other methods), and shed light onto the nature of the 
ground state structures in the neighborhood of the critical point. The main results 
obtained by this approach are confirmed by finite-size scaling analyses of numerical 
calculations based on the free-fermion method. 
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1. Introduction 

Quantum phase transitions differ from thermal transitions in that they happen at 
zero temperature, induced by quantum fluctuations ultimately inherent to all physical 
systems. As these transitions occur at zero temperature, they can be verified even in one- 
dimensional systems, which are especially sensitive to the effects of both quantum and 
geometric fluctuations, being thus ideal laboratories for the study of exotic phenomena 
associated with the presence of disorder. 

Arguably the simplest spin model to undergo a quantum phase transition is the 
quantum Ising chain, described by the Hamiltonian 



in which the a x ' z are Pauli matrices, Jj > represents a ferromagnetic bond between 
the longitudinal components of the spins at sites j and j + 1, and hj denotes a transverse 
field acting on the spin at site j. In the uniform limit (Jj = J, hj = h > 0), there is a 
quantum phase transition at h = J, between a ferromagnetic phase dominated by the 
nearest-neighbor bonds J, and a paramagnetic phase, dominated by the fields h. In this 
case, the model corresponds to the extreme anisotropy limit of the two-dimensional Ising 
model. Indeed, the quantum transition at h — J belongs to the Onsager universality 
class, whose paradigm is precisely the thermal transition of the two-dimensional Ising 
model. 

In the presence of disorder represented by randomly chosen couplings Jj and hj , the 
model in Eq. [T] exhibits surprising behavior. By extending a strong-disorder real-space 
renormalization group (SDRG) scheme introduced by Ma, Dasgupta, and Hu jTJ [2] to 
study the random Heisenberg chain, Fisher [3J H] was able to show that the nature of the 
quantum transition is fundamentally changed with respect to that of the uniform limit, 
being now characterized by a brutal distinction between average and typical behavior 
of various physical quantities, such as correlation functions. At criticality, there is a 
modification of the dynamic relation between length and time scales, which assumes an 
activated form (d la Arrhenius), in contrast with the usual power law. Sufficiently close 
to criticality, there are Griffiths phases [5], in which some thermodynamic properties 
are still singular. Fisher also provided exact calculations of a series of other properties 
for the random system, many of which are not known in the uniform limit, such as the 
scaling form for the magnetization as a function of an ordering field at the critical point. 

Partially similar effects are verified in the presence of aperiodic but deterministic 
couplings. This kind of aperiodicity is inspired by analogy with quasicrystals [6], 
which exhibit symmetries forbidden by traditional crystallography, corresponding to 
projections of high-dimensional Bravais lattices on low-dimensional subspaces. In one 
spatial dimension, aperiodic couplings can be defined from sequences of letters obtained 
from substitution or inflation rules, such as the one associated with the Fibonacci 
sequence, a — > ab, b — > a. Iteration of this rule yields a sequence abaababa . . . with no 




(1) 
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characteristic period. By associating different letters with different bond values J a and 
J;, (or with different field values h a and hb) one obtains an aperiodic quantum Ising chain. 
Any aperiodic sequence induces geometric fluctuations gauged by a wandering exponent 
(j, which describes the growth of these fluctuations as the chain length increases. A value 
of oj = | would mimic the geometric fluctuations induced by random couplings. 

For aperiodic quantum Ising chains, there are exact results obtained by RG 
treatments, unrelated to the SDRG scheme, and valid at the critical point El E]- 
These results allow the inference of the relation between length and time scales, and 
show that aperiodic sequences characterized by a wandering exponent u> = \ induce, 
at the critical point, effects analogous to those produced by randomness. Additionally, 
they confirm the validity of a heuristic criterion formulated by Luck |10j in order to 
gauge the effects of aperiodicity on ferromagnetic phase transitions. According to this 
criterion, aperiodicity is relevant, i.e. it is capable of changing the critical behavior of 
the quantum Ising chain, whenever it is associated with a wandering exponent u > 0. 
The marginal case u = may lead to nonuniversal critical behavior, dependent on the 
coupling values. 

In the present paper, we adapt the SDRG scheme to study the critical behavior of 
quantum Ising chains with aperiodic couplings. We show that this adaptation allows 
us to investigate the system both at the critical point and in its neighborhood, and to 
obtain presumably exact analytical results for a series of critical exponents, such as those 
associated with the bulk magnetization, the critical correlations, and the correlation 
length, which are not easily accessible via other methods. This allows a comparison 
with numerical results obtained for some of these exponents by Igloi et al. from a 
mapping of the quantum Ising chain to a free-fermion problem [TT| |12j. We note that 
essentially the same SDRG approach has been employed to study the scaling form of the 
entanglement entropy in aperiodic quantum Ising chains, although only at the critical 
point plE]. 

This paper is organized as follows. In Section [2] we review the SDRG approach 
introduced to deal with the random quantum Ising chain, as well as the associated low- 
energy physics. In Section [3] we introduce the kind of aperiodic structures in which 
we focus on this paper, as well as the notion of geometric fluctuations, discussing their 
connection with the Harris-Luck criterion for the relevance of aperiodicity to the critical 
behavior of the quantum Ising chain. Section [4] presents our extension of the SDRG 
approach to deal with aperiodic quantum Ising chains with couplings chosen from a 
family of aperiodic sequences whose geometric fluctuations can be tuned from irrelevant 
to relevant. In Section [5] we discuss the behavior of quantum Ising chains with couplings 
following Rudin-Shapiro sequences, whose geometric fluctuations mimic those associated 
with random couplings. The final section presents our concluding remarks. 
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2. Random quantum Ising chains and the strong-disorder RG approach 

The random quantum Ising chain corresponds to the extremely anisotropic limit of 
the McCoy- Wu model [T5l [16], and its very interesting low-temperature properties 
have been elucidated by Fisher (31 0], who argued that an adaptation of the strong- 
disorder renormalization group (SDRG) approach of Ma, Dasgupta and Hu [U E] yields 
asymptotically exact results for various physical quantities. (For a review of the SDRG 
method, as applied to the random quantum Ising chain and to many other problems, 
see Ref. [T7J.) 

The basic idea of the SDRG scheme is the iterative elimination of high-energy 
degrees of freedom, associated with locally strong couplings. For sufficiently strong 
randomness, in which there are very few strong couplings and many weak couplings, 
this is easily defined. In a finite chain, there is always a strongest coupling. If this 
coupling is a field, say h 2 , most likely the neighboring bonds J\ and J 2 are much 
weaker, and the spin a 2 will tend to align with the transverse field, giving negligible 
contributions to magnetic properties. If we are interested in the low-energy behavior, 
the spin o~ 2 , along with the field h 2 and the bonds J\ and J 2 , can be eliminated from 
the system. Nevertheless, virtual excitations will induce an effective coupling between 
the x components of the neighboring spins <j\ and 03, given, up to second order in 
perturbation theory, by 

If the strongest coupling is a bond, say J 2 , the spins o~i and o~ 2 that it connects are most 
likely under the action of much weaker fields h\ and h 2 , and thus will tend to respond, 
at low energies, as a single effective spin, under the action of an effective transverse field 
created by virtual excitations and given, again up to second order, by 
7 h x h 2 

hi = -j— ■ (3) 

■J2 

The bond J 2 and the fields hi and h 2 , along with the spins o\ and a 2 , are then replaced 
by the effective spin, with an effective magnetic moment, acted on by the field h\. 

In any case, the overall energy scale is reduced, since both effective couplings h\ 
and J\ are always weaker than any of the original couplings. As the transformations 
described above are iterated, this leads to a renormalization of the probability 
distributions of both bonds and fields, as well as lengths and magnetic moments. The 
fixed-point distributions can be found in closed form [H ITS] , and from these many low- 
energy physical properties can be derived. The nature of the ground state is strongly 
modified by the presence of disorder. While there is still a well-defined critical point, 
which corresponds to Pfeuty's exact result |19] . 



In hi = In Ji, (4) 

the critical relation between length (I) and energy (Q) scales departs from the power 
law valid in the uniform chain, having now the activated form 

fi-expf-v^), (5) 
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where Iq is some nonuniversal length scale. Most strikingly, in the neighborhood of the 
critical point the model is still gapless and the dynamic scaling behavior follows 

n ~ r*w, (6) 

in which z (5) is a dynamical exponent dependent on the distance to the critical point, 
measured by 

5 ~ In Ji — In hi. (7) 

This gapless noncritical regime arises from the existence of rare but arbitrarily large 
regions of spins which, due to statistical fluctuations, are locally in the "wrong" phase 
as compared to the whole chain, yielding the so-called Griffiths singularities |5]. 

Interestingly, the above results seem to be asymptotically valid not only for strong 
disorder, but for any amount of randomness present in the system. The reason is most 
easily seen at criticality, where the probability distributions of effective couplings become 
more and more singular around the origin, rendering the perturbative expressions for 
the effective couplings essentially exact. In the following sections, we show that also in 
the presence of relevant aperiodic modulations it is possible to obtain asymptotically 
exact results, although the nature of the ground state turns out to be quite distinct 
from that of the random case. 



3. Aperiodic sequences and the Harris-Luck criterion 



Substitution rules defining aperiodic sequences have the general form 

{ai ->■ Wi , i 6 {1,2, . . . ,n}, 

in which Wi is a "word" formed by letters in the set {a{\. The sequence is built by 
iteratively applying the rule to an initial letter. For instance, the celebrated Fibonacci 
sequence is obtained by applying the rule 

a — > ah 
h — > a 

to an initial letter a. 

Various properties of the resulting infinite sequence (obtained after an infinite 
number of iterations of the rule) are related to the substitution matrix 

/ fll fl2 " ' fin \ 
fl\ fll ■ ■ ■ fin 

\ fnl fn2 ' ' ' fnn ) 

in which f\j denotes the number of letters a« in the word Wj. For instance, the length 
of the sequence after k iterations of the substitution rule is asymptotically given by 
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Ci > being the largest eigenvalue of M. Also, the ith entry of the corresponding 
eigenvector (normalized so as the sum of all entries equals unity) gives the relative 
frequency di of the different letters in the infinite sequence. 

For the purposes of the present paper, the most relevant property derived from 
the substitution matrix are the fluctuations of the relative frequency of a letter, after k 
iterations of the rule, with respect to that of the infinite sequence. These fluctuations 
are explicitly defined as 



(*) 
91 



N® - diN k 



where N® denotes the number of letters a,- after k iterations of the substitution rule. 

,0 



It can be shown that g k , for all i, is asymptotically governed by the second largest 
eigenvalue (in absolute value) of M, 



in which 



to 



2 



InCi 

is the so-called wandering exponent of the aperiodic sequence. A negative value of 
to indicates that fluctuations are bounded, and become negligible as the length of 
the sequence becomes larger and larger. On the other hand, a positive u> leads to 
unbounded fluctuations, as in the case of a sequence built from independently and 
randomly distributed letters, for which u — |. Finally, the marginal case of oj = 
usually leads to fluctuations which grow logarithmically with the sequence length. 

When couplings in a ferromagnetic model are chosen in a nonperiodic fashion, the 
Harris-Luck heuristic criterion [TUt [20] makes use of the concept of the fluctuations of 
those couplings in a length scale defined by the correlation length to point out under 
which conditions aperiodicity changes the critical behavior with respect to that of the 
uniform system. In the case of a quantum Ising chain with bonds chosen from some 
aperiodic sequence, such fluctuations are proportional to the g® defined above, and 
the criterion states that aperiodicity will be relevant, i.e. the critical behavior will be 
changed, whenever the wandering exponent satisfies 

cu > 1 - i = 0, (8) 

in which v — 1 is the correlation-length critical exponent of the Onsager universality 
class. 

In the following section, we analyze the effects of bonds chosen from a family of 
aperiodic sequences with wandering exponents ranging from u) = — 1 to u — > 1. 

4. A family of aperiodic sequences 

Consider the family of aperiodic sequences generated by the substitution rule [21] 

a ^ ahk , b k = bb---b (9) 
b -> a ^ ' 
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(a) 



(b) 



(c) 



Figure 1: Renormalization-group transformation as applied to the Fibonacci chain. The 
original chain is depicted in (a), with weak J a bonds indicated by dashed (blue) lines 
and strong Jb bonds by solid (red) lines; spins sit on the open circles, under the action 
of the original transverse field h. After the spins connected by strong bonds are joined, 
the chain takes the form shown in (b), with open circles now denoting effective fields h 
and solid circles representing original fields h. Decimation of the spins under the action 
of h leads to the sequence of effective couplings shown in (c), which again corresponds 
to a Fibonacci chain. 



with k a positive integer, so that k — 1 corresponds to the Fibonacci sequence, illustrated 
in Fig. [l](a). The substitution matrix has eigenvalues 

Cf = \ ± \ViT±k, 

and thus the wandering exponent, 
ln |C | In A; 

u k = + = — p - 1, 10 

is a monotonically increasing function of k, with U\ = —1, U2 = 0, W3 = 0.317. . ., etc. 
According to the Harris-Luck criterion, we thus expect that aperiodicity is relevant for 
k > 3, with k = 2 representing a marginal case. 

The sequences generated by Eq. [9] contain clusters of k letters b, separated by 
clusters of 1 or k + 1 letters a. In order to implement a strong-disorder RG treatment 
of the corresponding quantum Ising chain, we have to calculate the effective bonds and 
fields associated with those clusters. Although we start with a uniform field, the RG 
scheme will generate different effective fields for different clusters, and this will be taken 
into account in the discussion below. 

We first look at the case of n spins Uj (i — 1, 2, . . . , n) under the action of a field h , 
forming with spins o~i and a r a cluster connected by n + 1 bonds J <^ ho] see Fig. |2](a). 
Furthermore, suppose that the fields hi and h r acting on o~\ and o~ r are small compared 
to /io- Physically, one expects that the n spins in the cluster tend to line up with the 
field ho, and thus can be decimated out of the system. Applying perturbation theory up 
to order n + 1 leads to the conclusion that fluctuations of the cluster induce an effective 
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Figure 2: Renormalization of couplings and decimation of spins in the aperiodic quantum 
Ising chain, (a) Elimination of spins connected by n neighboring strong fields h . (b) 
Formation of a effective spin from a cluster of n spins connected by strong bonds Jq. 



coupling 



T71+1 

J=-T- 



between o~i and o~ r . 

A more delicate situation corresponds to a cluster of n spins connected by strong 
bonds Jo, while the field h acting on each of the n spins, as well as the bonds J\ and J r to 
the rest of the chain, are much smaller; see Fig. [2^b). However, this can be easily dealt 
with by resorting to a duality transformation [22], which interchanges the roles of bonds 
and fields. In the transformed space, the situation is precisely the one of the previous 
paragraph, so that the result can be immediately written down. At low energies, the n 
spins behave as a single magnetic moment, under the action of an effective field 
r h n 

h=~^, (12) 

with an effective magnetic moment 

\i = nfi, (13) 

H being the magnetic moment of each of the n spins. 

Getting back to the sequences generated by Eq. |9j suppose we look at intermediate 
values of the field, so that J a h J^. (The case J a 3> leads to a renormalized 
chain in which the roles of weak and strong bonds are interchanged, as long as the 
system is sufficiently close to criticality. Thus, there is no loss of generality in assuming 



<h ^> Ja-) At energy scales Q such that h < O < J&, Eqs. 12 and 13 imply that each 
cluster of k+ 1 spins connected by bonds Jb behaves as a single magnetic moment, under 
the action of an effective field 



i 



h=^r, (14) 
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with an effective magnetic moment 

jt = (k + l)fi. (15) 

As h < h, and J a < h, the largest coupling is now provided by the original field h, as 
illustrated in Fig. [ljb) for the case of the Fibonacci chain. Further reducing the energy 
scale, so that J a < Q < h, the clusters of k + 1 bonds J a , containing k spins under the 
action of the field h, can be decimated out, leaving effective bonds 

J" +1 



Jn 



(16) 



as shown in Fig. [jjc). All original fields h are decimated in this step, so that all active 
spins at this energy scale are under the action of the same field h. Now the largest 
couplings in the chain are provided by the remaining clusters of isolated J a bonds, 
which we can interpret as effective bonds 



Jb Jai 



(17) 



since with this choice the sequence of effective bonds exactly reproduces the original 



sequence. Equations 14 through 17 define an RG step. 



In general, there is no guarantee that the effective couplings satisfy the condition 
J a < h < J b , which is required to iterate the RG procedure. However, we shall assume 
that the condition is satisfied, and explore the consequences. We first introduce a change 
in notation, rewriting the RG equations as 



kC>+i) 



■J a 



k+1 



J 



Cj+i) 



° n 1 



k+l 



J, 



(i+i) 



Li) 



(18) 



(19) 



in which Ja , jjf , and /j,^ represent the values of the parameters after j iterations 
of the RG procedure, with J a ^ = J a , = J b , = h, and jj^ = 1. Next, we define 
the ratios 



J, 



U) 



and 



■J a 



U) 



(20) 



which both remain smaller than unity as long as the condition Ja^ < < J ; 



r(i) 



is 



fullfilled. From Eqs. [181 and 19 we obtain recursion relations for the ratios, 



[r®]' 



and 



fe+1 



which can be written in matrix form by taking the logarithm of the previous equations, 
yielding 



l nr (i+ 1 ) 
InsO+i) 




(21) 
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Defining 



\ v j) 



In r-O') 
In 



and 




(22) 



and expressing the matrix as Tk = UfeD^U^ 1 , in which the diagonal matrix consists 
of the eigenvalues of H\ and the columns of the matrix contain the corresponding 

as 



DfcU," 1 \v s ) = 



right eigenvectors, we can rewrite Eq. [2T 

\ u j+i) = u fc 1 \ v j+i) 
from which we readily obtain 

\ Uj ) = K \u ) => \ Vj ) = UfcD^ 1 \v ) . (23) 

The eigenvalues of the matrix Tk are related to the eigenvalues of the substitution 
matrix by 



A; 



±\2 



while the matrix Ut. can be written as 



1 




Using the above results, we can invert and expand Eq. 23 to obtain 



and 



In 



lns (i) 



1 



{KY*k + {KY^ 



(KY C* k - (\ k y Ck 



(24) 



(25) 



with 



a± = C±lnrW + lns (0) . 
We thus see that the behavior of r^' and s^) for j ^> 1 is in general dominated 



> 0. If a k , which enters the coefficient of 



by the largest eigenvalue of A^ > A A 
(X^Y in both lnr^ and lns^, is positive, we asymptotically have lnr^ <C and 
lns^ ^> 0, and thus the effective fields are much smaller than the effective bonds. This 
corresponds to the ferromagnetic phase. Conversely, if a k is negative, the asymptotic 
behavior corresponds to lnr^ 3> and lns^ <C 0, the effective fields are much larger 
than the effective bonds, and the system is in the paramagnetic phase. The critical 
point must then correspond to the condition = 0, which, in terms of the original 
parameters, leads to a critical field 



h 



f-< k j 



crit 



jd a jd-b 



(26) 



in which d a and db = 1 — d a can be shown to be the fractions of letters a and b in 
the infinite sequence produced by Eq. |9j This last result agrees with the exact result 
derived by Pfeuty [TU] for the critical condition of the quantum Ising chain. That the 
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approximate RG equations used here are able to lead to the exact result in Eq. [26 



must be related to the fact that the recurrence relations 14, 16 and 17 preserve the dual 
symmetry of the quantum Ising chain. 

It is also useful to determine the effective lengths of the various couplings as the 
RG scheme is iterated. We follow Fisher [4] and assume that fields and bonds share 
lengths evenly, so that the initial bonds have lengths 
has length 



"ft 



From Eqs. 



and 



19 



(o) _ g(p) _ i an j ^ e j n jtj a i field 



and by looking at Fig. [IJ we see that these 



lengths satisfy recursion relations which can be written in matrix form as 



V 



-b 
"ft 



V 






k 



k 


k + 1 



\ 



( 6P \ 



J 



Li) 
b 

U) 
ft 



( & \ 



\ w ) 



I k + 1 
1 


Thus, all effective lengths grow asymptotically with the largest eigenvalue of the matrix 
Sfc. It turns out that the eigenvalues of S fc are 1, A^T and A^ , so that 

*W~^~^~(AJ)''. (27) 

In the following subsections, we analyze the consequences of the previous 
expressions for the behavior of the system both at the critical point and in its 
neighborhood. 



4-1. Properties at the critical point 



The recursion relations 18 and 19 can be iterated indefinitely at the critical point. As 



then a k = 0, the behavior of the ratios and s^' is determined only by the smallest 
eigenvalue Ar of the matrix T k , 



In r 



a, 



and 



In s 



(KY 



(28) 



Cfc ~ (k (k ~ (k 

Precisely at the critical point, there are exact real-space RG calculations [7J El [9] 
that provide information on the dynamic scaling between energy levels and length scales. 
These can be used to gauge the reliability of the SDRG approach developed here. 



4-1.1. Irrelevant aperiodicity. For the Fibonacci chain, which corresponds to k — 1, we 
have A^ ~ 0.382, so that both and asymptotically approach unity. This means 
that the perturbative recursion relations become poorer and poorer approximations as 
the iterations proceed, suggesting at the same time, in agreement with known results 
(see [7] and references therein), that the critical behavior should correspond to that of 
the uniform chain. In other words, Fibonacci modulations are irrelevant for the critical 
behavior of the quantum Ising chain, and the relation between energy and length scales 
at criticality follows the usual dynamical scaling form 

with a dynamic exponent z — 1. 
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4-1.2. Marginal aperiodicity. For k = 2, we have \ k = 1, so that both and 
remain constant along the iterations, being expressed in terms of the original bonds as 

so that the ratio between the effective bonds is constant, 

jU) j 

(i) _ ^_ _ (i) (j) _ _ 

This corresponds to a line of critical points, associated with different values of the ratio 
p = J a / Jb, and gives rise to a nonuniversal dynamic exponent z, as already known for 
other aperiodic sequences [7JIE1IH], an d for quantum Ising chains with correlated disorder 
[23J. Indeed, if we use the values of the strong bonds at each step of the RG procedure 
to estimate the gaps IX,- between different energy levels, we get 

n, ~ J b U) = Jt l) = $^Jt 1] = pJt 1] = fa, 



while Eq. 27 leads to a characteristic length 



Expressing j in terms of £j, we then obtain the dynamic scaling form 

n 3 ^q z(p \ (29) 

with a dynamic exponent 

/ \ In P , % 

= ( 3 °) 

As in the case of the quantum spin-^ XX chain with marginal aperiodic bonds 
\Z4\ [23] . the dynamic exponent coming out of the SDRG approach corresponds to the 
leading term in the p 1 expansion of an exact result which can be obtained by other 
methods [HI El EE]- m the present case, the exact result turns out to be 

Z e, (P) = g ' ( 31 ) 

which can be obtained from the corresponding result for the quantum Ising chain with 
bonds chosen from the period-doubling sequence generated by the substitution rule 

A ->■ AB 
B -> AA ' 

Distributing couplings and ^> Ja according to the above sequence, and applying 
an intermediate field h, the SDRG approach leads to a sequence of bonds J a and J& 
given by the rule in Eq. [9] with k = 2. At criticality, these bonds are related to the 
original ones via 

5/3 



/ 

T — A T — T 

2/3' 

B 
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so that 



T / T \ 2 / 3 



Thus, the exact result for the period-doubling chain |8j|9] translates into Eq. 31 which 



for p <C 1, clearly reduces to Eq. 30 



4-1.3. Relevant aperiodicity . For k > 3, the wandering exponent Wfc is positive. 
According to the Harris-Luck criterion, this corresponds to relevant aperiodicity. 

It is easy to check that, under the condition k > 3, is larger than unity. Since 



we start with r^ < 1 and < 1, we have < 0, and Eq. 28 implies that both 
and s^') become arbitrarily close to zero as the RG scheme proceeds. This means 
that the perturbative recursion relations become asymptotically exact, as does the RG 
treatment itself. This is analogous to what happens for the random quantum Ising chain 
[31 H]. It is then not surprising that the activated dynamic scaling form 

^■~exp[-(V01> ( 32 ) 
with £ p ~ |lnp|~ 1//aJ a length scale that depends on the original couplings, is exactly 



reproduced by the SDRG treatment, as it can be readily checked from Eqs. 18, 19, 26 



and [27} Thus, the dynamic exponent is formally infinite. This result is analogous to 
what is verified for the critical spin-^ XX Z chain with relevant aperiodic bonds [231 [25]. 

Notice that u here plays the role of the "tunneling" exponent ip in the random 
quantum Ising chain [27J, which owes its name to the need of high-energy virtual 
excitations in order to flip spin clusters. 

4-2. The neighborhood of the critical point 

In order to study the system in the neighborhood of the critical point, we set 

h = h CTit (1-5), \5\ < 1, (33) 



with 5 measuring the distance to criticality. Taking into account Eq. [26j and the fact 

and 7 are 



33 



that now In J, = In h CT it, one can see that the definitions of 5 given by Eqs. 
equivalent for h ~ fo C rit- 

It is clear from the previous discussion that, for 5^0, the RG process has to be 
interrupted whenever the condition 

max{lnr (i) ,lns (j) } = (34) 

is reached, which happens for a value of j we denote by jg. By using Eqs. [24[ and [25| 
together with the initial values 

r<°> = \ = p d ° (1 - 6) and = ^ = p db (1 - 5)" 1 , 

Jb h 

in which we again make use of the ratio p = J a / Jt between the original bond values, we 
can show that for \5\ «C 1 the condition in Eq. [34] leads to 



. In | ,5 1 Hi5| . . 

1 ^^KK) +c ^^^KK-y (35) 
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c being a constant which is different in the paramagnetic (5 < 0) and in the ferromagnetic 
(5 > 0) phases. 

This last result allows us to obtain analytical expressions for critical exponents, 
as discussed below. Of course, the reported results are not valid for k = 1, where the 
critical behavior corresponds to that of the uniform chain. 



4-2.1. The paramagnetic phase. The paramagnetic phase is characterized by the fact 
that after jg iterations all effective bonds are smaller than the effective fields. Thus, the 
length of the largest effective bonds at this point of the RG process provides an estimate 
of the correlation length. 



Substituting Eq. [35] into Eq. [27] we obtain 

e-^- to)*~i*r, 

with a correlation-length exponent 



lllA;h (36) 



In A^T — In A fc 



As A^ = (Cfc") 2 ; this last result can be rewritten as 

1 



(37) 



1 -u k 

We thus see that the correlation-length exponent at the aperiodic critical point saturates 
the Harris-Luck criterion in the form 

v > — !— , (38) 
l — UJ 

for all sequences in the family generated by the substitution rule in Eq. |9j In particular, 
for the marginal case of k = 2 we recover the exact result v = 1 EJ E]. We note that 



Eq. 38 is analogous to the Chayes bound for the correlation length critical exponent in 
the presence of disorder |28j . 

We can also study the behavior of the energy gap as a function of 5. From the 



discussion in Subsection |4.1[ we see that the energy gap vanishes at the critical point, 
since Qj — > as j — > oo; see Eqs. 29 and 32 On the other hand, in the paramagnetic 



phase the value of the effective field for j = jg, being much larger than any other 
effective coupling, yields an estimate of the energy gap between the ground state and 
the first excited state of the infinite chain. We thus have 

n(S) ~ ~ \5\ z{p)v , 
for marginal aperiodicity {k = 2), and, 

n (S) ~ fcfo) ~ exp (-C IC^) > (39) 
for relevant aperiodicity (C being a nonuniversal constant). The above expressions can 



be derived from the recursion relations in Eqs. [18] and [19] with the help of Eq. 35 
Notice that Q (5) — » as 5 — > 0, but the noncritical gap is always finite, implying the 
absence of Griffiths singularities in the presence of aperiodic fluctuations generated by 
substitution rules [T2] I2T)]. 
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4-2.2. The ferromagnetic phase. In the ferromagnetic phase, after j$ iterations of the 
RG transformation the effective field becomes smaller than both effective bonds. An 
estimate of the average spontaneous magnetization is provided by the average magnetic 
moment of all the spins still active at this point of the RG procedure (i.e. those that 
have not been decimated under the action of strong local fields in the previous RG 
steps). Since the sequence of bonds remains invariant up to this point, the estimate is 
simply given by 

m x ~ n j6 ^ s) — , 

in which n, is the fraction of active spins after j iterations of the RG transformation. 



From Eq. [19j we see that 

n& = (k + iy^ , 

yielding 



Substituting Eq. 35 then leads to 
with a critical exponent 

Notice that, if we follow Fisher [31 H] and define at the critical point an exponent <f> 
relating the effective magnetic moment to the energy scale, 



where f2 p is a nonuniversal constant, we find 

p = v (l - 44) , 

just as for the random quantum Ising chain 
4-3. Critical correlation functions 

Now we determine the asymptotic behavior of the pair correlation function associated 
with the order parameter, 



C xx (0 = <« + />, 

in which the bar denotes average over all lattice sites i. At criticality, this correlation 
function follows a power-law scaling form, 

C xx (I) ~ ^. 
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In the uniform limit, and thus also for the Fibonacci chain, the exponent rj takes the 
value | [29]. Below we calculate r] for k > 2. 

The elements associated with the strongest bonds jff^ in each RG step are clusters 
containing k + 1 spins. Each of these clusters, if decoupled from the rest of the chain, 
has a doubly degenerate ground state (in the absence of a field), corresponding to all 
spins pointing along the +x or the — x direction. If we denote these states by |\? ), it 
is clear that 

(tt* laf]^) = ±1, 

in which i labels the spins inside a given cluster. 

Owing to the inflation symmetry associated with couplings generated by the 
substitution rule in Eq. [9j these clusters are either grouped into superclusters of k + 1 
clusters or are separated from each other by k + 1 weak bonds J a , as illustrated for 
k = 2 in Fig. [3]^a). The superclusters, on the other hand, are either grouped into 
"supersuperclusters" containing k+1 superclusters or isolated from each other by clusters 
and weak bonds, and so on. Thus, there appears a hierarchy of (k + l)-spin clusters, 
yielding a self-similar structure of the correlation function, and leading to the appearance 
of characteristic distances where the correlation is particularly strong. 

In order to analyze this structure, we recourse to the simplest approximation, 
in which we assume that spins appearing in the same cluster at some step of the 
RG scheme are maximally correlated, while the remaining correlations are negligible. 
Looking at Fig. |3](a), we see that, in the limit of a very large chain (where boundary 
effects can be ignored), the average correlation between spins separated by the distance 
l = £^ + = 1 , corresponding to the real distance between nearest-neighbor spins 
connected by a coupling Jb, is given by 

1 n' 0> m <0) 

in which N is the total number of spins, is the number of Jb bonds (riV b being 
the number of clusters with k + 1 spins), and the multiplicative k factor comes from the 
fact that in each cluster there are k pairs of spins separated by the distance Iq. In fact, 
we have 



C xx (l' Q ) 



k + 1- g ^ N b 



(o) 



x 



k N ' 

since in each cluster there are k + 1 — go pairs of spins separated by the distance l' Q = golo, 
go € {1,2, ... ,k}. However, the strongest correlation corresponds to the value g = 1, 
coming from spins separated by the distance 1' = Iq. 

From Fig. |3j(b), we see that an estimate of the average correlation between spins 



separated the distance l\ = £^ + = k + 1 is given by 



C xx (h) 



b ~ <-h 

N 



1 N b V 



x (k + 1) = (k + 1) 
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Figure 3: Estimating correlation functions for the case k = 2. 
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with denoting the number of effective bonds Jj and the multiplicative factor k + 1 
coming from the fact that each effective spin represents k + 1 real spins in this step 
of the hierarchy. Notice that inside each cluster there are in fact contributions to the 
average correlation between real spins separated by distances l[ = g\l\ ± goto, where 
gi G {1, 2, . . . , k} and go G {0, 1, . . . , k}. (However, these correlations are strongest for 
l[ = li.) Since l\ — kl — 1, l\ — (k — 1) Iq — 2, . . ., l\ — Iq — k, this gives corrections 
to the average correlations C xx (Iq), C xx (21q), . . ., C xx (Mo) between spins separated by 
short distances in the real chain. 

Figure [3tc) leads us to conclude that the average correlation between real spins 

h 



separated by the distance l 2 = £^ + I? — k 2 + 3k + 1 is approximately given by 



C xx (h 



1 



AT(2) 

x(A; + l) 2 = (A; + l) 2 ^ r . 



since now every effective spin represents [k + l) 2 real spins. In fact, there are 
contributions to the average correlations between real spins separated by distances 
1' 2 = g 2 h ± ± 9oh, in which g 2 G {1, 2, . . . , k} and g x , g Q G {0, 1, . . . , k}. (Again, the 
strongest contribution comes from spins separated by the distance 1' 2 = Z 2 -) Here there 
are corrections to correlations between spins separated by distances of the order of li, 
but not of the order of l Q , since l 2 — kli — kl = l\ corresponds to the minimum value of 
l> 2 . 

Generalizing the above discussion, there are strong correlations associated with 
distances around lj = £^ + £^ ~ {^t) 1 > w ifh 

c xx (i J ) = (k + iy^ r , 

Njf being the number of effective bonds in a chain with ^> 1 sites. Actually, 
there are contributions to correlations between real spins separated by distances 

l 'j = 9jlj ± 9j-llj-l ± • • • ± fi-o^o, 

with gj G {1,2,..., A;} and g^j G {0,1,..., k}, so that the strongest contribution 
appears for E = lj, and the minimum value of E is of the order of lj-i- Thus, there are 
no corrections for correlations between spins separated by distances of the order of l^ 2 
and below. 

Owing to the invariance of the sequence of effective bonds along the RG steps, the 
number of effective J b bonds scales as 

Af ~ f . 

It is then clear that we can estimate the long-distance scaling behavior of the average 
strong correlations as 

c-^-^t^-i (41) 
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with an exponent 

In v 

This result is in contrast with the one obtained for the random quantum Ising chain, 
for which critical correlations decay with distance as C xx (I) ~ l~ 2 P/ v ^ due to the 
probabilistic nature of cluster formation. 

We note that the prediction in Eq. 41 is valid only for the characteristic distances lj. 
Correlations between spins separated by generic distances (distinct from the E defined 
above) are expected to be much weaker. In fact, as shown explicitly for the XX chain 
in Ref. [25J (Appendix B), these weaker correlations are proportional to the effective 
couplings, and thus follow essentially the same scaling form as in Eq. 29 (32 ) for marginal 
(relevant) aperiodicity. 



4-4- Numerical calculations 

The analytical results derived in the previous subsections can be checked by numerical 
calculations based on the well-known mapping of the quantum Ising chain onto a free- 
fermion system \30\ [31]. Here we use the recipe described in Ref. |32j . The numerical 
work then involves diagonalizing matrices of order 2L for chains containing L spins. As 
here we work with open boundary conditions, this corresponds to L — 1 bonds. 

When dealing with finite chains, we have to take into account that substitution rules 
produce aperiodic systems containing many different subsequences of the same length 
L — 1. We then have to perform a weighted average over all O (L) such subsequences, 
with weights given by the frequencies with which a given subsequence appears in the 
infinite sequence. For arbitrary L, in order to determine these frequencies we iterate 
the substitution rule until a very large sequence of size C (with C ranging from ~ 10 4 
to ~ 10 6 ) is obtained. We then calculate the finite-size estimates of the corresponding 
frequencies, and extrapolate to infinite C using standard algorithms |33j . In order to 
minimize fluctuations due to the discrete scale invariance exhibited by substitutional 
sequences, the number of bonds, L — 1, is always chosen to coincide with the length of 
the aperiodic sequences after an integer number of iterations of the substitution rule. 

We work here in the so-called canonical ensemble of disorder, according to which 
we study the critical properties of the system by performing numerical calculations for 
a value of the transverse field corresponding to the critical field of the infinite sequence. 
The other possibility, which defines the microcanonical ensemble, would correspond to 
setting the transverse field to values calculated from the densities of letters a and b in 
each particular subsequence of finite length, which are usually distinct from the infinite- 
sequence value. It has been shown that the canonical ensemble yields the same critical 
exponents as exact and RG calculations |32] . 

The finite-size behavior of the bulk critical magnetization is described by 

m x (L) ~ L-?/\ 
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Figure 4: Rescaled magnetization profiles for k = 3 at the critical point, for different 
chain lengths L. The values of the exponents and v were set to their SDRG predictions, 
yielding /3/u ~ 0.169. Calculations were performed for fixed ends, (erf) = (af) = +1, 
and the curves result from averaging over all distinct subsequences of length L, followed 
by mirror symmetrization around the position I = L/2. 



for systems of size L. Within the free-fermion approach, this scaling behavior can be 
checked by calculating magnetization profiles (i.e. local expectation values of the spin 
operators erf, / = 1, 2, . . . , L) in finite chains with prescribed boundary conditions [32J. 
Here we focus on fixed boundary conditions, corresponding to h\ = and hi = 0, and 
look at the resulting profiles when both end spins are in the state \af) = |+). 

For k = 3, a case in which aperiodic modulations are relevant, we expect the SDRG 
prediction for f3/v to be valid for any coupling ratio p = J a / Of course for J a j Jf, close 
to unity there will be a large crossover length L x , below which the system will appear 
uniform, and the SDRG result should be recovered only for L > L x , while for smaller 
values of p the SDRG should be evident even for relatively small chain lengths. This is 
confirmed by the free-fermion calculations, as shown in Fig. |l]for J a / Jb — 1/10. Notice 
that the results for finite L (L between 41 and 1160) collapse onto two different master 
curves, depending on whether the subsequence size was produced after an odd or an even 
number of iterations of the substitution rule. This is related to the fact that the deviation 
of the number of letters a (or equivalently b) with respect to the expectation value in 
the infinite sequence oscillates in sign between consecutive iterations of the substitution 
rule (the eigenvalue £2 of the substitution matrix is negative). Furthermore, it should be 
noted that none of the master curves follows the conformal invariance prediction |3"4"t 133] 

m x (l/L) = (o-?)~[Lsm(nl/L)]-P /u , 

as expected from the strongly anisotropic scaling along the space and time directions 
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Figure 5: Critical bulk magnetization for k = 2 (actually, for the related period-doubling 



chain described in 4.1.2) as a function of the chain length L. The upper dashed curves 
correspond to the Onsager behavior L~^' v ', with (3/u = 1/8, while the lower dashed 
curve indicates the SDRG prediction, (3/v ~ 0.2075. 



implied by Eq. 32 



In the marginal case, k = 2, we only expect the SDRG predictions to be satisfied 
for coupling ratios approaching zero. Indeed, as shown in Fig. [5j a plot of the finite- 
size behavior of the bulk magnetization at criticality (estimated from the corresponding 
value in the middle of the chain) shows a clear change from the Onsager universality 
class to the strong-modulation behavior as the coupling ratio is reduced. 

The behavior of the critical pair correlation function C xx (/), as obtained from the 
free-fermion method, is illustrated in the main panel of Fig. [6] for k = 2 and a small 
coupling ratio J a /Jb = 1/100. Notice the existence of characteristic distances for which 
spins are strongly correlated, due to the presence of strong effective bonds, and of a 



ultrametric structure produced by the cluster hierarchy discussed in Subsection 4J3 
The inset shows plots of the typical correlation, estimated by 



C»(/) = exp (\n(afa- 



+i 



which exhibits nonuniversal behavior and a much smoother dependence on the spin 
separation I, confirming that, as in the random quantum Ising chain, there is a 
decoupling between average and typical behavior. As contributions from the rare 
strongly correlated spin pairs are washed out by the logarithm in the calculation of 
Ctyp, typical correlations behave just as the weak correlations discussed in the end of 
subsection 14.31 
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Figure 6: Main panel: Average pair correlation function versus spin separation for k = 2. 
The red dashed curve shows the Onsager behavior, while the black dashed curve indicates 
the SDRG prediction, which is valid at the characteristic distances where correlations 
are strong. Inset: Typical correlation function for various values of the coupling ratio. 
In all cases, numerical calculations correspond to chains with 1366 sites. 



5. Rudin-Shapiro sequences 



5.1. The four-letter sequence 

The family of sequences analyzed in the previous section has the property that, at the 
criticality, there is only one type of spin cluster connected by strong bonds at a given step 
of the RG process. A more complex situation is produced by choosing bonds according 
to the four-letter Rudin-Shapiro (RS) sequence, generated by the substitution rule 

a — > ab 
b — > ac 
c -)■ db 
d — > dc 



(43) 



The associated substitution matrix has eigenvalues 0, ±a/2, and 2, so that the wandering 
exponent is uj = |, the same produced by uncorrelated randomness. In the resulting 
infinite sequence all letters appear with the same frequency, and no clusters of repeated 
letters are present. 

We choose couplings such that J& = pj a , J c = p 2 J a ,and Jd = p 3 Jd, with < p < 1. 
Based on Pfeuty's result |19j . we anticipate a critical field 

3 

^crit = P^Jai 
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Figure 7: Renormalization-group transformation as applied to the four-letter Rudin- 
Shapiro chain. The original chain is shown in (a): J a bonds are represented by thick 
solid lines (red); by solid lines (yellow); J c , by thick dashed lines (blue); J^, by 
thin dashed lines (cyan); open circles mark the position of the spins, under a uniform 
field h ~ h CTit . Spins connected by J a or bonds form clusters, while spins connected 
only by J c or Jd are decimated. These processes give rise to the effective chain in (b), 
characterized by 6 different values of both bonds and fields, respectively indicated by the 
numbers below and above the chain. Upon further lowering of the energy scale, effective 
spins under the grey boxes cluster together, while crossed circles denote effective spins 
decimated by the action of strong effective fields. As shown in (c), apart from a minor 
boundary effect, the next effective chain preserves the structure in (b). 
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Jl 


J CI 


hi 


= h 2 /J b , 


/'i 


= 2, 


■h 


= Jdl 


h 2 


= h 2 /J a , 


/'2 


= 2, 


J 3 


= J c Jd/h, 




= h 3 / J a Jb, 


/'3 


= 3, 




= J 2 c J d /h\ 


hi 


= h 4 / J a J b -, 


fl 4 


= 4, 


■h 


= JcJ 2 Jh\ 


h 5 


= hyj 2 a j b , 


/'5 


= 4, 


J 6 


= JtJl/h\ 


h 6 


= hS/JtJl 


/'6 


= 5. 



so that, in order to look at the system in the neighborhood of the critical point, we 
focus on the condition Jd ^ J c ^ h ^ J b ^ J a . At energy scales Q such that 
Jd < Jc < ^ < h < Jb < J a, all spin pairs connected by J a or J b bonds behave as single 
magnetic moments, whereas spins connected to their neighbors at both sides by only J c 
or Jd bonds are effectively removed from the chain by the action of the transverse field; 
see Fig. [71(a). The resulting effective chain is composed of 6 different bonds and fields, 
whose values (in decreasing order of intensity), and those of the corresponding effective 
magnetic moments, are given by 



(44) 



As illustrated in Fig. [7]^b), weak bonds appear next to strong fields and vice- versa. 
Writing again 

h = h cl . it (l-5), (45) 

which defines 5 as a measure of the distance to criticality, it can be seen that for \5\ <C 1 
and p <C 1 the couplings are ordered according to 

J 1 ~ h x > J 2 ~ h 2 > J 3 ~ h 3 > J 4 ~ h A > J 5 ~ h 5 > J 6 ~ h 6 . 

Sufficiently close to the critical point, upon further reduction of the energy scale, this 
sequence of bonds and fields is mapped onto itself by the bond recursion relations 



4 j) /h?, 
J® J®/ [h? 
J^J^/hfh^ 
jU) jjJ) rO') / 





-1) 


- J U 

— J A 


Jt 


-1) 


~ J U 




-1) 


~ J U 




-1) 


~ J U 




-1) 


~ J U 


J U ~ 


-1) 


~ J U 



(46) 



J. 



(j) 



J, 



Jk JJ / h\ 

2 



(3) 



6 



Itii 



h 



(?) 



the field recursion relations corresponding to the dual versions (hi -H- Jj) of the above 
equations, thus defining the RG transformation. We also write recursion relations for 
the effective magnetic moments, 



u ij+1) - u U) 

/4 i+1) = riP + rfP + n?, 

/4 i+1) = 



(47) 
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as well as for the lengths of the effective bonds and fields, 



,0+1) 
/i+i) 
/i+i) 

1 J A 

/?0+l) 
l J6 



?0) 
'J4> 

$-0) i 

J2 ^ 

?0) i 

'J2 ^ 

s+ 



./2 



+ 



48 + 

?0) i 

-J3 ^ 

?0) , 

-J4 ^ 

»0) 



?0) 

-hi ' 

+ 2 ^(i) 



\/3 



+ 



-J5 
»0) 



hi ' 
0) , 
hi 



>0) 

h2 ' 



(48) 



?0) i tfO) 



J5 



+ 



(i) +2^ + 



-J3 



'J4 



-J6 



J6 

0) , 9 /i) 

hi ^ ^ c h2 ) 



>Cj) 

h2 ' 



with analogous expressions obtained by interchanging the labels J and ft. Here, lj) or 
Zj^ represent the length of the effective bond Jj or field ftj, i G {1,2,..., 6}, at the jth 
iteration of the RG transformation. 

The recursion relations in Eqs. [46] and their dual counterparts are valid as long 



as hi > J2 and > h^K The first condition eventually ceases to be valid in 
the ferromagnetic phase, while the second one fails in the paramagnetic phase. To be 
concrete, let us focus on the ferromagnetic phase. From the dual symmetry of the 
recursion relations (under the interchange of bonds and fields), similar results are valid 
in the paramagnetic phase. 

As in the previous section, it is convenient to introduce a coupling ratio 



,0) 



J. 



U) 



h\ 



(J)' 



(49) 



whose value allows us to check the validity of the condition hf > J^ 1 ■ In order to 
obtain a closed set of recursion relations, it is also necessary to define the additional 
ratios 



rii) 



,0) 



J, 



0) 



h 



(J) 



and 



,00 



h) 



J, 



or 



ze{l,2,...,6}. 



Substitution of Eqs. 46 into the above expressions, and definition of the ratio 



t 



U) 



,0')„0) 

'1 s 4 



,(i+i) 



(50) 



(51) 



lead to a set of recursion relations which in matrix form can be written as 
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2 J 
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v j+1 ) = T \vj) , 
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in a notation analogous to that of Eqs. [21] and [22} From Eqs. [44pT| we see that the 
initial value of the vector \v) takes the form 



\ v o) 



/ 21n(l- 


S) 




21n(l - 


<)) 




41n(l - 


<)) 




61n(l - 


<)) 




61n(l - 


5) 




81n(l - 


5) 




\ 2 In p - 8 In 


(1- 


5) / 



The eigenvalues of the matrix T are Ai = 4, A2 = 2, A3 = A4 = 1, and 
A5 = A6 = A7 = 0. Denoting the respective right eigenvectors by \4>i), « € {1,2,..., 7}, 
\v ) can be expanded as 



\v ) = -8 In (1-6) \<p 1 )+2\np\ 



^}--ln(l 



5) + N> 



where \wq) belongs to the kernel of T. Thus, from \vj) = TP \vq) and from explicit 
expressions for \4>i), \<p 2 ) and | ^3) , we obtain, for i G {1,2, . . . 6}, 

lnr? = (frXi + gtXi) In (1-5), 
IntCi) = -8Ailn(l-5) + 2A|lnp, 

with the constants and all positive. It is then clear that at the critical point (5 = 0) 
we have r 



Li) 



1, i.e., h 



(j) 



j\ at all RG steps, whereas asymptotically approaches 
zero. This indicates that the RG scheme becomes asymptotically exact, in agreement 
with the relevant nature of aperiodic modulations characterized by a wandering exponent 
co = I > 0. Accordingly, the dynamic scaling can be shown to assume the activated 



form 



exp 



with Ej a characteristic length scale and l p a function of p only. This is exactly the form 
obeyed by the dynamic scaling behavior of the critical random quantum Ising chain 

0131. 

In the ferromagnetic phase (5 > 0), we have < 1, as expected, but the RG 
scheme is no longer valid once (or equivalently t^; see Eq. 51) reaches unity. This 
happens for 

~41n(l-<5)~ 



lnt (r) = 



hi 



3 



In p 



In 5 



In 



ln(A 2 /A! 



From this step on we need to redefine the RG transformation, since now 
It turns out that after one step of this new RG transformation all effective fields are 
smaller than all effective bonds. Thus, critical properties can be estimated, in the same 
way as in the previous section, by looking at effective quantities after js ~ j* iterations 
of the RG scheme. 



(52) 



,(i) 
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Figure 8: Rescaled average magnetization profiles for the four-letter Rudin-Shapiro 
sequence, with a coupling ratio p = jq- 



The matrix connecting effective lengths in consecutive RG steps, built from Eqs. [48 
and their dual counterparts, has as its largest eigenvalue Ai = 4, so that the correlation 
length can be estimated as 

with a correlation length exponent 

In Ax 1 

^ = = 2 = . 

In Ai — In A2 1 — OJ 

again saturating the Harris-Luck criterion as written in Eq. |38} 



From the recursion relations of the effective magnetic moments, Eq. [47J we similarly 
extract a characteristic magnetic moment scaling as 

in which $1 = 3.24698 ... is the largest eigenvalue of the matrix connecting j/^^ j to 
' ^ Thus, the spontaneous (bulk) magnetization should scale as 



/V 



t*36 



#1 

Ai 



.IS 



6? 



with a critical exponent 

ln(VA0 



/3 



0.300902.... 



ln(A 2 /A!) 

In order to check the validity of the above results, we again resort to exact 
diagonalization of finite chains, based on the free-fermion method. As shown in Fig. 
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[8j these calculations are fully compatible with the above prediction for the ratio (3 jv. 
Notice that now there is a clear collapse of the curves only for the bulk magnetization, 
and the presence of four rather than two distinct bond leads to the disappearance of the 
second master curve visible in Fig. |4| 

5.2. The two-letter sequence 

The results derived in the previous subsection are valid as long as the four bonds assume 
distinct values. However, by making the identifications = J a and = J c we obtain 
a binary RS sequence, and the properties of the corresponding quantum Ising chain can 
in principle be different. For instance, numerical calculations and scaling considerations 
based on a mapping to a directed walk led Igloi et al. [12J to propose for this binary 
sequence a correlation-length critical exponent v = 4/3 and a ratio /3/u = 0.160 (5). 

For the quantum Ising chain with couplings following the binary RS sequence, with 
Ja — pJc, < p < 1, the critical field is given by 

^crit = P 2 Jc- 

When applying the SDRG scheme, we thus focus on the condition J a <C h <C J c . As 
shown in Fig. [9]^a), the original chain has clusters with 1 to 4 strong bonds. At energy 
scales Q such that J a < h < Q < J c , this gives rise to 4 different effective magnetic 
moments {p n } and fields {h n }, with 

h n+1 

fi n = (n + l)fJi and h n = — n 6 {1,2,3,4}. (53) 

At this point, all bonds are equal to J a , and the largest energy scale in the chain is 
provided by the original fields h. Thus, at energy scales such that J a < Q < h, all spins 
under the action of h can be decimated, inducing 4 different effective bonds 

jn 

J " = /^I' {1,2,3,4}. (54) 

The resulting effective chain is depicted in Fig. |9^b), and corresponds to a sequence of 
4 different bonds and fields. 

Exactly at the critical point, the bonds and fields are ordered such that 

J 1 = hi > J 2 = h 2 > J 3 = h 3 > J 4 = h A , 

with 

n + l 

Jn = P 2 Jc- 

Notice from Fig. ^b) that now, contrary to what happens for the four-letter RS chain, 
there are clusters of 1 to 3 strong bonds J\ or fields h\, and in some of these clusters both 
strong bonds and strong fields are present. Of course, this picture remains valid in the 
neighborhood of the critical point, making it harder to write recursion relations, since 
the perturbative treatment becomes more involved. However, we can take advantage of 
a mapping between the critical quantum Ising chain and the quantum XX chain [36J 
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Figure 9: Renormalization-group transformation as applied to the two- letter Rudin- 
Shapiro chain. The original chain is shown in (a): J a bonds are represented by dashed 
lines (blue), while stronger J c bonds appear as solid lines (red); open circles mark the 
position of the spins, under a uniform field h ~ /i C rit- Spins connected by J c bonds 
form clusters, while spins connected only by J a are decimated. These processes give 
rise to the effective chain in (b), characterized by 4 different values of both bonds and 
fields, respectively indicated by the numbers below and above the chain. Upon further 
lowering of the energy scale, effective spins under the grey boxes cluster together, while 
crossed circles denote effective spins decimated by the action of strong effective fields. 
As shown in (c) the next effective chain preserves the structure in (b). 
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and of the results derived for the aperiodic quantum XX chain (see Ref. |25| ) to write 
the bond recursion relations 



■ J i 



Cj+i) 



J 



-Cj+1) 

2 

•O'+l) 
3 

O'+i) 



J, 



2 



U) 



j(J) rii) tU) jU) I uU) 

■ J i J 2 J 3 J i I 'n 



(55) 



(j) jU) 



J 9 t^ f . 



(?) 



/ 



with dual recursion relations for the effective fields, the function f (J,h) being explicitly 
given by 



f(J,h) = (j 2 + h 2 )-^. 

The above recursion relations define an RG transformation which maps the effective 
chain onto itself as long as 

> h { i ] and hf > 4 j) . (56) 

Recursion relations for the lengths of the effective couplings take the form 
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with the initial values 

At the critical point, / h/j 
recursion relations in Eqs. 55 recover the multiplicative structure satisfied by the 
previously discussed aperiodic chains. Then, we can show by the method described 
for the four-letter RS sequence that again the SDRG approach becomes asymptotically 
exact, predicting a dynamic scaling form 



cxp 



consistent with the value u = \ of the wandering exponent of the binary RS sequence. 

The analytical treatment of the recursion relations in the off-critical regime is 
hindered by the presence of the function f(J,h). It is possible to perform an 
expansion on the distance 5 to criticality in order to determine at which step of the 



RG transformation the condition in Eq. 56 ceases to be valid. However, not only is the 
resulting equation impossible to solve in closed form, but also the new RG transformation 
which we have to impose does not seem to lead to an invariant sequence of effective 
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Figure 10: Dependence of the lengths Is (proportional to the correlation length) and lj* 
on the distance to criticality S, for a coupling ratio J a /Jc=\- 



couplings. Thus, we have to recourse to a numerical implementation of the SDRG 
scheme. 

In order to reach the largest possible system sizes, we start at a given distance 
to criticality 5 (defined as in Eq. 45) and apply the recursion relations 55 (and their 
dual counterparts) up to the point at which min I J*f \ 



< max 



keeping 

track of the effective lengths via Eqs. [571 Let us denote the effective couplings at 



this point by 



r(i*) 
•Jn 



and < hn 



A j \h 



Since up to this point the sequence of effective 



couplings is invariant under the RG transformation, it can be generated by applying 
the rules leading to Eqs. 53 and [54] to a chain with bonds following the binary RS 



sequence. To the resulting effective couplings, we then assign the values 

tin 



T U*) 
On 



and 



and iterate numerically the SDRG scheme by sweeping through the lattice and 
progressively eliminating the high-energy degrees of freedom. For a binary RS sequence 
with 2 20 ~ 10 6 letters, the corresponding number of { J n } and {h n } couplings is around 
5 x 10 5 , and correcting for the lengths of the effective couplings leads to a real chain 
containing 6 x 10 12 spins for 5 = 10 -5 . 

We implement this numerical procedure, in the paramagnetic phase, up to the point 
at which all bonds become smaller than all fields. In Fig. [lO we plot, for p — |, the 
scaling behavior of the length is of the strongest bond at that point, as a function of 



at 



the distance to criticality. Also shown is the length £j* of the strongest bond j{ J ' 
the point beyond which the sequence of effective lengths is no longer invariant. We see 
that is nicely follows a power law £$ ~ \o~\~ 2 , and thus we conclude that, using ig as an 
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Figure 11: Rescaled average magnetization profiles for the two- letter Rudin-Shapiro 
quantum Ising chain. Main panel: rescaling with the exponents calculated by the SDRG 
approach. Inset: rescaling with the numerical estimate of Ref. 1121 . The dashed line in 



the main panel corresponds to the conformal-invariance prediction, Eq. 58 



estimate of the correlation length, we obtain 

with an exponent v — 2. 

This result is in contrast to that obtained by Igloi et al. from a scaling argument for 
the surface magnetization of finite chains, based on an equivalence of the problem with a 
directed walk. As already mentioned, these authors predict a correlation-length critical 
exponent v — g. They also provide an estimate of the ratio fijv ~ 0.160, by adjusting 
numerical results for critical magnetization profiles of small chains (containing between 
9 and 65 spins), with a ratio between weak and strong bonds given by J a /J c = 1/16. 
Incidentally, they find the surprising result that the data for such small chains can be 
well fitted by the predictions of conformal invariance, although the system is clearly not 
conformally invariant, in view of the strongly anisotropic scaling along the space and 
time directions. 

We then decided to use the free-fermion method to look at larger chain sizes, 
and smaller bond ratios, in order to check whether we could fit the numerical results 
with our prediction for /3/v, assuming, on the grounds of universality, that the ratio 
(3/v — 0.150451 . . . found for the four-letter RS chain would still be valid for the two- 
letter RS chain. Figure 11 shows the results of free-fermion calculations of critical 
magnetization profiles for chains containing between 65 and 513 spins, and J a / J c = 1/50. 
Despite the small difference between our prediction for 0/u and that of Ref. |12| . the 
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Figure 12: Finite-size scaling plot of the lowest gap calculated numerically for off-critical 
(paramagnetic) quantum Ising chains with couplings following the binary Rudin-Shapiro 
sequence. Parameters correspond to u = \ and v — 2. The dashed lines are proportional 
to (\5\ L X I V Y" V . Values of \b\ range from 1(T 4 (1(T 3 for the inset) to 8 x KT 1 (8 x 10° 
for the inset) for all system sizes. 



inset shows that our prediction yields a better data collapse. Although the results shown 



in Fig. [TT] are definitely incompatible with the conformal-invariance form 

m x (l/L) ~ [L sin (irl/L)]~ p/v , (58) 

we do observe a reasonable agreement with such form for small chains (not shown), as 
in Ref. |12| . This agreement becomes poorer as the coupling ratio is reduced. 



Furthermore, the results in Eqs. 32 and 39 can be used to check directly the 



prediction v — 2, from free-fermion calculations of the lowest gap Ql {8) for chains with 
L spins, at a distance 5 from criticality, via the finite-size scaling ansatz 

-lnQ L (5) = f (\8\ L 1 ^) , 

in which / (x) is a scaling function predicted to behave as 

/(x)~x— 

for x ^> 1, corresponding to large system sizes at a fixed 5, and to approach a constant 



as x — > 0, which corresponds to \S\ — >■ 0. As seen in Fig. 12, we obtain a very good data 
collapse for binary Rudin-Shapiro chains with coupling ratio J a /J c — 1/10, and lengths 
ranging from L—l = 64 to L—l = 4096, by setting oj = | and using the SDRG prediction 
v — 2. (The calculations were performed in the paramagnetic phase, corresponding to 
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5 < 0, for which the lowest gap is obtained from the smallest positive fermion energy. 
In the ferromagnetic phase this vanishes exponentially with the system size, and the 
lowest gap is obtained instead from the second smallest positive fermion energy.) The 
main plot in Fig. 12 shows the rescaled lowest gap produced by the subsequence of 
L — 1 bonds (L = 65, 257, 1025, 4097) obtained, starting from a single letter a, by 
iterating the substitution rule in Eq. 43, and imposing Jt = J a and = J c . For the 
sake of illustration, the inset shows the corresponding results obtained by averaging the 
logarithm of the lowest gaps for all subsequences of L — 1 bonds cut from an infinite two- 
letter Rudin-Shapiro sequence. (Due to computational costs, chains with L = 4097 were 
not available for averaging.) Results for the average lowest gap, as well as for the average 
second gap, give qualitatively the same scaling picture (not shown), in agreeement with 
SDRG predictions. 



6. Discussion and conclusions 



In this paper, we have presented an extension of the strong- disorder renormalization 
group approach for the study of the low-energy critical behavior of the quantum Ising 
chain with couplings chosen from aperiodic but deterministic sequences, generated by 
substitution rules. 

For the cases of marginal and relevant aperiodicity, in which the critical behavior 
is expected to depart from the Onsager universality class, we have been able to obtain 
analytical results for various critical exponents, such as those related to the correlation 
length, the bulk magnetization, and the spin-spin correlation functions at criticality. 
Interestingly, we always find that the correlation-length critical exponent is given by 
v = (1 — w) -1 , a result which barely satisfies the Harris-Luck criterion |TD S 20J around the 
aperiodic fixed point, which takes the form v > (1 — u)^ 1 . We have also derived results 
for the off-critical lowest energy gaps, which reflect the absence of Griffiths singularities 
for the present class of aperiodic systems, as expected from previous arguments [T2| 126]. 
All our results have been fully confirmed by numerical calculations on finite chains, 
based on the free-fermion methods. 

The basic physical picture offered by the SDRG approach, derived for a family 
of two-letter substitution rules, has been shown to remain valid for more complicated 
sequences, represented by the two- and four-letter Rudin-Shapiro sequences. According 
to this picture, sufficiently close to the critical point, as the energy scale is progressively 
reduced, the effective chains generated along the RG process exhibit a self-similar 
structure, in which the sequence of effective couplings remains invariant and their values 
decrease, whereas the magnetic moments of the effective spins and the distances between 
those spins grow exponentially. This picture holds up to a certain step of the RG process, 
beyond which the ground state of the chain settles into either the ferromagnetic or the 
paramagnetic phases, characterized respectively by the dominance of bonds or fields. 

The self-similar nature of the effective chains is helpful in understanding the absence 
of Griffiths singularities in the systems studied here. At all length scales, self-similarity 
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guarantees that geometric fluctuations are entirely governed by the wandering exponent 
u, setting a maximum size for locally critical clusters, as already noted in j26]. This 
is in contrast with the uncorrelated random quantum Ising chain, where the wandering 
exponent only gauges the growth of average geometric fluctuations, there being a nonzero 
probability for the occurrence of arbitrarily large clusters of spins which are locally in a 
different phase than the bulk. 

Finally, the results reported here, concerning critical exponents and the nature 
of the ground-state phases, remain valid for aperiodic XY chains, described by the 
Hamiltonian 

H = ~\ E J < K 1 + 7i) «+i + (1 - 7.) - £ 

i 

as long as the anisotropy parameters 7j, < |7j| < 1, are all positive or all negative. 
Under these conditions, the SDRG scheme is straightforwardly adapted, and it can be 
shown |37j that the anisotropics j$ flow towards | -y^ | = 1 as the RG scheme is iterated, 
so that the low-energy behavior corresponds to that of the quantum Ising chain. 
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